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1. Introduction 

Recent data reported by the PAMELA 1 and ATIC 2 experiments, as well 
as cosmic ray lepton results from FERMI 3 and H.E.S.S. 4 collaborations 
sparked a series of papers explaining these results either by contributions 
of electron positron pairs to the local interstellar cosmic ray (CR) spectrum 
by dark matter (DM) or pulsars. Focusing here on pulsars, we argue that at 
the present, our knowledge about particle acceleration at pulsars as well as 
of the local Galactic CR propagation is still limited, i.e. the recent results for 
CR electrons and positrons constrain pulsar and propagation models. We 
will thus not present another attempt to explain the data by contributions 
of pulsars to the local CR lepton flux but rather to highlight the caveats in 
doing so. 



2. Pulsars as Particle Accelerators 

Modelling the acceleration of particles at pulsars and subsequent injection of 
these particles into the interstellar medium, one may look at three scenarios: 
young pulsars (after the breakup of their pulsar wind nebulae), 5 mature 
pulsars 6 and millisecond pulsars. 7 
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2.1. Energy Budget 

Before discussing the acceleration of particles by pulsars, it is instructive to 
look for the available energy budget of a single pulsar. The energy reservoir 
of a pulsars is its rotational energy l/2f2 2 , where I is the moment of inertia. 
The spin down power L$d = of a pulsar is given by 8 

W*) =L sdfi {l + t/r )-^ , (1) 

with the characteristic decay time r = Po/((n — l)Po) and n = 3 for a 
magnetic dipolc field of the pulsar. 

With Pq, Pq being the period and its derivative at pulsar birth and 
assuming the magnetic field of the pulsar does not decay, the spin down 
power at birth can be written as 



L s d,0 — 



4tt 2 IPP 



pi 



(2) 



It is easily seen from Eq. 2 that, given a small enough birth period, one can 
get an arbitrary large energy reservoir. 10 ms can be seen as a reasonable 
lower limit for P , while initial periods in the range «50-150ms are not 
uncommon. 9 It is estimated that up to 40% of the pulsars may be born 
with periods in the range 100-500 ms. 10 



2.2. Source Spectra 

Several models have been put forward to describe the acceleration of CR 
leptons by pulsars (see e.g 11 and references therein). 

The simplest approach is to assume that some fraction / part of the spin- 
down power Lsr>(t) is transferred to particles, which follow a power-law 
spectrum with index a and exponential cut-off. The acceleration of leptons 
at 7-ray pulsars has been modeled by a number of authors. Whereas 12 
assume that the positron spectrum follows the 7-ray spectrum, it is assumed 
by 6 that mature pulsars with ages larger than 100 kyr inject their mono- 
cncrgctic wind into the interstellar medium. 

From Eq. 1 it is clear that the available spin-down power (and thus 
the available energy that can be transferred to particles) is largest when 
the pulsar is still young. In fact, integrating Eq.l over time for n = 3, 
it becomes clear that at t = r the pulsar dissipated half of its available 
rotational energy. One can therefore expect that the majority of particle 
acceleration at a given pulsar takes place for t > To . What is indeed observed 
are nebula of relativistic particles around many young pulsars, i.e. pulsar 
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wind nebulae (PWN). There is also evidence, that these highly relativistic 
particles are confined in the PWN. 13 

Assuming the particles from the pulsar are reaccelerated at the pulsar 
wind shock and then contained in the PWN until its breakup, 5 model the 
particle spectrum injected into the interstellar medium by the two middle 
aged, nearby pulsars B0656+14 and Geminga. We generalize this model by 
assuming a broken power law particle spectrum at the pulsar shock 14 

f k'{t)(E/E B y 2 for E>E B 

Q'(E,t) = { , . (3) 

[ k'^iE/EB)- 1 for E<E B 

keeping E B a free parameter. The spectrum at the shock is normalised by 
the available energy deposited into particles 

Q'(E,t)EdE = / part L SD (t). (4) 
The fraction of the spin down power deposited in particles 



J E, 



/part =/?/(! + CJ) (5) 

can be expressed in terms of the magnetisation parameter 

a{t) = 0.003 (i/lkyr) 3/2 . (6) 
With e = ri/r s hock = 0.001 . . . 0.1, the maximum particle energy at shock 



is 14 



£ max = een v / a/(a + 1)a/L sd /c. (7) 

k = 3 is the compression ratio at the shock and e the elementary charge. 

The evolution of the particle spectrum in the PWN is governed by syn- 
chrotron losses 

- Q'(E,t) = A (B PWN (tfE>Q(E,t)) (8) 

where the mean PWN magnetic field is assumed to decay with time 

SpwnM = 1200/(1 + t/kyr) 2 \pG]. (9) 

Under these assumptions, the particle spectrum in the PWN at time of 
breakup T is: 

Q{E,T)= / Q'(E ,to)EoE- 2 <d(E -E niin )e(E niax -Eo)dt ,(10) 
Jo 

with 

Eo=E/(e £° B? WN {t r ) 2 dt' + 1^ 
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3. Propagation of Cosmic Ray leptons 

The propagation of CR leptons in the Galaxy is dominated by diffusion and 
energy losses due to synchrotron radiation in the Galactic magnetic field. 
The differential density TV satisfies the equation 

f)N f) 

— -S = V.(kVN)- — (b E*N), (If) 

where bo is determined by the Galactic magnetic field. 
In the case of a pulsar, the source function is given by 

S = Q(E,t)5(f-f SOUIce ), (12) 

where Q(E, t) is given by a model for particle acceleration at pulsars. Above 
« 4 GV, the diffusion coefficient has an energy dependence of the form 

k = k {E/l GeV) a , (13) 

where a — 0.3 ... 0.7 and ko are determined by fitting observed secondary 
to primary and radioactive secondary data. Depending on the model, ko is 
in the range 0.006kpc 2 Myr~ 1 -0.2kpc 2 Myr" 1 (e.g. 15 ). 

Although obtaining ko in this way is a well established procedure, we 
want to remark the following. First, fitting secondary to primary (e.g B/C) 
data only yields a correlation between k and the halo height H. 16 This 
degeneracy can be broken by looking at radioactive to stable secondary ra- 
tios (e.g. 10 Be/Bc). In doing so, one compares life times due to radioactive 
decay with the times scale of escape, the latter depends on the diffusion per- 
pendicular to the Galactic disk. Thus the fitting of the observed chemical 
composition of CR by a propagation model measures primarily the propaga- 
tion perpendicular to the Galactic disk, which not necessarily has to be the 
same as that in the Galactic plane. Second, to the best of our knowledge, all 
studies determining ko by fitting the observed chemical composition so far 
use ID or 2D models. Nevertheless, it has been shown, 17 that the primary 
component of locally observed CR flux may change at least by a factor of 
two depending on the local SN history, given these objects are indeed the 
main sources of Galactic CR. The flux of the secondary component on the 
other hand shows little or no variation. Thus, the observed secondary to 
primary ratios will depend on the local SN history. It may be necessary to 
focus on tertiary to secondary ratios when working with leaky box, ID or 
2D propagation models. 
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4. Modelling the Contribution of Pulsars to the local 
Cosmic Ray Lepton Spectrum 

If the source function of electrons and positrons accelerated at a given 
pulsar and also the distance to this object is known, one can calculate the 
contribution of these particles to the locally observed CR spectrum. As 
discussed above, the magnitude of the diffusion coefficient in the Galactic 
plane is not well determined. Also, the spectrum of the particles injected 
into the interstellar medium strongly depends on the model used. 

Focusing on particle acceleration at young pulsars as, according to Eq. 1, 
these objects promise the largest output of high energy leptons, we follow 
the model outlined in the previous section. According to pulsar physics, 
we have to allow for a large range in parameter space. The pulsar birth 
period can be in the range w 10 ms to 500 ms, Eb, is most probably in 
the range 10-1000 GeV, and the PWN life time T smaller 100 kyr. Within 
this parameter space in addition to the uncertainties in k , it is possible to 
produce a wide range of particle spectra at earth which can fit the data. 
In fact, looking at Geminga or B0656+14 individually, the available CR 
lepton data imposes constraints to the parameter space of our model. 

5. Experimental Discrimination between Dark Matter and 
Pulsar origin 

If indeed the observed increase in the positron fraction is due to the pulsars 
B0656 and Geminga, the fact that those two pulsars are located about 
7 degree apart on the sky in the outer Galaxy, we expect an anisotropy 
towards the galactic anticenter, as discussed in. 5 However, for DM we expect 
either an isotropic sky, or an anisotropy towards the galactic center region. 

In the diffusion model, the anisotropy in the CR flux is given by 18 
6 = 3k \VN\/(cN). For instantaneous injection at time t% and distance r» 
for a diffusion coefficient as given by Eq. 13, one arrives at 

S = A ri bo ( Q _ i) e* (V- 1 - (E/ ((t - U) b E + l))"- 1 ) ~* , (14) 

5 find an anisotropy of a few percent above 10 GeV. We remark that one 
expects for E < 10 GeV a solar cycle dependent anisotropy that is due to 
modulation of CR electrons and positrons. 

6. Summary 

The ATIC report of an excess just below 1 TeV resulted in a wide inter- 
est in the interpretation of the signal. This signal is also supported by the 
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PAMELA detection of an increase in the positron to electron ratio towards 
the energies where ATIC saw their excess above the general cosmic ray 
lepton (electron and positron) spectral background. Interestingly, this ef- 
fect was not seen by the detector on the FERMI satellite, 3 nevertheless 
the latest FERMI results indicate a flattening of the CR electron spectrum 
above « 50 GeV ', hinting at an additional leptonic component with a hard 
spectrum. 3 If we assume that the ATIC peak is real, then we are left with 
two competing explanations: (i) Kaluza-Klcin type particle, or, (ii) a lo- 
cal source(s) of electrons and the best candidate for the latter is a pulsar 
(with its PWN) since both the Kaluza-Klcin and the pulsar interpretations 
require both electrons and positrons, probably in equal numbers. 

Whereas the Kaluza-Klein explanation required the introduction of a 
boost factor of 200 in the volume production rate to explain the ATIC 
signal, 2 the pulsar/PWN interpretation also involves a number of free pa- 
rameters, such as the birth period and interstellar diffusion coefficient be- 
tween the source and earth. At the moment it seems as if we are left with too 
many free parameters and uncertainties to make an unambiguous statement 
about the origin of the ATIC (or at least the PAMELA) signal. More work 
needs to be done on DM candidate particles, as well as the multiwavelength 
nature of pulsars and their wind nebulae. 
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